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Abstract: We show that the 0{A) ambiguity in the pole mass can be fixed in a natm'al 
way by introducing a modified nonperturbative V-scheme momentum space coupfing ay (q) 
where the confining contributions have been subtracted out. The method used is in the 
spirit of the infrared finite coupfing approach to power corrections, and gives a nonpertur- 
bative definition of the 'potential subtracted' mass. The short distance expansion of the 
static potential is derived, taking into account an hypothetical short distance linear term. 
The magnitude of the standard OPE contributions are estimated in quenched QCD, based 
on results of Liischer and Weisz. It is observed that the expansion is not yet reliable at the 
shortest distances presently measured on the lattice. 
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1. Introduction 

Historically, the pole mass M and the heavy quark potential V{r) were among the first 
quantities where renormalons Q have been discussed in a physical context in QCD. Latter, 
the connection of the 0(A) ambiguity in the pole mass [||, ^ with a corresponding ambiguity 
in the coordinate space potential Q was pointed out. It was observed js], ^] that the 
leading renormalon contribution cancels in the total static energy Estatic = 2M + V{r), a 
physical quantity which should be free of ambiguities. This cancellation is a non-trivial 
finding. Indeed, one might have expected that the pole mass and the static potential 
should be separately well defined: for instance, in the Schrodinger equation, the quark 
mass normalizes the kinetic energy. Furthermore, although the potential appears to be 
nonperturbatively defined only up to an arbitrary constant (in particular only the force is 
the quantity free of ambiguity in lattice calculations), it is difficult to maintain the view 
that the arbitrary normalization of V{r) implies an arbitrary normalization of M, which 
nevertheless would follow from the non-ambiguity of the static energy if there were no 
independent way to fix the normalization of either the mass or the potential. In this paper 
I suggest that there is in fact a natural way to define unambiguously the pole mass at the 
nonperturbative level (at least as far as the leading renormalon ambiguity is concerned) even 
in a confining theory like QCD, by properly subtracting out the confining contributions to 
the self-energy, hence to fix also the 'constant term' in the potential. In Sec. 2, the definition 
of the 0{K) term in the pole mass is given, in term of a properly defined nonperturbative 
momentum space V-scheme coupling aY{q). The method used is in the spirit of the infrared 
(IR) finite coupling approach to power corrections [^. In Sec. 3, theoretical constraints on 
av{q) are reviewed. In Sec. 4 the short distance expansion of V{r) is derived, including the 
effect of an hypothetical linear short distance term, and the standard IR power corrections 
are estimated on theoretical ground. It is shown that present lattice data are not available 
at distances short enough for a reliable short distance analysis to be performed yet. 
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2. The nonperturbative pole mass 



To define the pole mass, one has to fix its weh-known renormalon ambiguity [0, |3[. I 
start from the result ^ that the leading IR contribution 5Mprp^j^ to the perturbative 
pole mass MpT (when expressed in term of a short distance mass like m = tuj^), is 
related (presumably to all orders of perturbation theory ||]) to the leading long distance 
contribution 6Vpt\ir to the perturbative coordinate space potential Vpt by the relation 



^MpT\jp{fif) = -^^VpT\IR{^J^f), (2.1) 



where 



^VpTimiH) = I 7^'^pt{q)- (2-2) 



VpT{q) is the momentum space perturbative potential, related to VpT{r) by Fourier trans- 
formation 

VpT{r) = I ^ exp{iq.r)VpT{q), (2.3) 

and /i/ is an IR factorization scale. Defining to all orders of perturbation theory a momen- 
tum space potential effective coupling av\pT{Q) by 

VpTiq) = -AttCf^^^^^P^, (2.4) 



eq.( |2.lD can be rewritten as 



Cp 

SMpT\ip{fif ) = / dq av\PT{q)- (2-5) 

^ Jo 



The right hand side of eq.( |2.5| ) is presumably ill-defined, since it involves an integration 
over the IR Landau singularity thought to be present in ay\pp{q), and represents (taking 
fj,f A) the 0{A) ambiguity in the pole mass. To solve this problem, one would be 
tempted, in analogy with the IR finite coupling approach to power corrections [^, to 
replace the perturbative effective coupling av\PT{l) inside the integral in eg. (p. 51) by the 
corresponding nonperturbative coupling ayiq) defined by 

V{q) ^ -AttCf^^, (2.6) 
q^ 

where this time V{q) is the Fourier transform of the full nonperturbative potential V{r): 



^(0 = / TTTTTs exp{iq.r)Viq). (2.7) 



d^q 
("2^ 

However, in a confining theory, V{q) either does not exist (e.g. if V{r) ~ i?logr -|- C for 
r — > oo), or is anyway too singular at small q (refiecting the singular large distance behavior 
of V{r)), making the integral in eq.(p.5|) (with the nonperturbative av{Q)) divergent at 
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g = 0. For instance, in the case of a linearly raising potential V{r) = 0{r) for r oo, 
one gets av{q) = 0{l/q^) for q —>■ 0. This observation suggests one should first subtract 
out the confining long-distance part of the potential to define a suitable nonperturbative 
coupling av{q)- To this end, the following procedure appears the most natural one: expand 
the potential around r = oo, and subtract from V{r) the first few leading terms in this 
expansion (including an eventual constant term) which do not vanish for r ^ oo. There is 
by construction only a finite number of such terms. Let us call their sum Vconf {'>')■ Then 
we have 

V{r) = Vconf{r) + 6V{r), (2.8) 

which, assuming the large r expansion can actually be performed, uniquely defines 6V{r), 
such that 5V{r) both for r ^ (from asymptotic freedom) and for r — > oo. It is clear 
that 5V{r) now admits a standard Fourier representation 

SV{r) = I expiiq.r) 5V{q), (2.9) 

and one can define the new nonperturbative coupling av{q) by 

5Viq) = -4ttCf^^. (2.10) 
q-^ 

One should note that the perturbative part of these quantities are preserved, namely 
5VpT{r) = VpT(r) and 6VpT{q) = VpT{q), since 5V differs from V by the Vconf term, 
which, viewed from short distances, appears as a finite sum of nonperturbative power-like 
corrections, invisible order by order in perturbation theory. Indeed, the terms occurring 
in perturbation theory should scale as 1/r, hence vanish for r oo, which excludes them 
from Vconf (r). Thus av\PT{q) = (^viPTiq) is the same as in eq.(|2.4|), i.e. ay and ay have 
identical perturbative expansions. 

As an example, consider the potential in quenched QCD (this is actually the only case 
where the analytic form of the r — > oo expansion is known in low orders). Theoretical 
expectations give the long distance expansion for r — > oo 

V{r)c:,Kr + C-~ + 0{\), (2.11) 
12 r 

Although the 0{l/r) term is not a rigorous result of QCD, since it has been derived within 
an effective bosonic string theory Q, it has been numerically confirmed ^ in high precision 
lattice simulations. We shall therefore assume that eq.(^TlD gives the correct large distance 
behavior of the static potential. It follows that 

Vconf{r) = Kr + C, (2.12) 

and one defines 

V{r) = Kr + C + 5V{r). (2.13) 
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In this case, the couphngs ayiq) (if it can be defined nonperturbatively, i.e. if C = as 
previously noted) and av{Q) just differ by a term, arising from the Fourier transform 
of the Kr piece. 

The prescription for the nonperturbative definition of the pole mass now reads as fol- 
lows. Introduce the 'potential subtracted' mass 

Mpsifif) = MpT - 6MpT\iR{jif), (2.14) 
and define the nonperturbative IR contribution to the pole mass by 

5MIR{^if) = -]^5VIR{^^f), (2.15) 

where 

5VIR{^^f) = [ 5V{q), (2.16) 

J\q\<tif 

which yields 

6MiRifif) = — / dq dviq), (2.17) 
Jo 

in complete analogy with eq. (|2.lD , (|2.2| ) and ( |2.5D . Then the pole mass is given by 

M = Mps{fif) + 6MiR{fif) + ..., (2.18) 

where the dots represent non- leading 0{l/m) IR contributions from higher order renor- 
malons, and the fif dependence approximatively cancels between the two terms on the right 
hand side. The interpretation of the prescription eq.( |2.18D is transparent: it says one should 
remove from Mrt its ambiguous IR part dMpj^^jR^Hf), as suggested in |5|, and substitute 
for it the corresponding nonperturbative (and non-ambiguous) IR contribution 5MiR{fif). 
One should note the similarity between eq. (|2.18| ) and the corresponding expressions in the 
IR finite coupling approach to power corrections [Q]. In the present context, however, the 
nonperturbative coupling is unambiguously identified. With the pole mass well-defined, 
the constant term C in the large distance expansion of the potential (eq.( p.H"|) ) is in turn 
fixed, since the corresponding constant term in the large distance expansion of Estatid^), 
which should be unambiguous and calculable, is 2M + C. 



3. Constraints on the nonperturbative av{q) 

Eq.(glTD and (|2l^) yield 6Vir) for r ^ oo, hence 5V{q) ~ -^^ for q ^ 0, 



which yields 

Cpdy(g = 0) = ^, (3.1) 

i.e. 6iv{q = 0) ~ 0.196, a rather small IR fixed point value. Substituting this value as a 
rough estimate of av{q) in the integrand of eq.( p.l7 ) gives 
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Cf 1 
6MjR{nf) ~ —av{q = ^) f^f = /^/' 



(3.2) 



which represents a correction of about 100 MeV for the range of quoted in Q for 
b-quarks. 

A more refined estimate is obtained by inputting the information about the 0{\) 



term in eg. ( 2.11 ), which was obtained in [|9[| from a fit to high precision large r lattice data 
and yields for r — > oo 



6V{r) 

with b ~ 0.04/m, hence for g — > 

6V{q) 



vr 1 IT b 
T2r ~ 12^ 



^1 
' 3 q2 



6 q 



and 



(3.3) 



(3.4) 



(3.5) 



Note that, since 6 > 0, av{q) increases from its IR value as q increases, hence must be 
non- monotonous in the IR region, since asymptotic freedom implies it should ultimately 
decrease to at large q. At /// = 1.2 GeV, the second term in the parenthesis in eq.(p^) 
represents a correction of about 40% to the IR value. Substituting eq.( |3.5| ) in the integrand 
of eq.( p.l7D gives 



which yields 6Mir{iij) ~ 120Mey for /x/ = 1.2 GeV. 



(3.6) 



4. Short distance expansion of the heavy quark potential 

In this section I show that, barring constant terms, the short distance expansion of the 



heavy quark potential can be obtained directly^ from eq.(2^), despite the singular behavior 



of V{q) at small q. Introducing again the factorization scale /ij, eq.(2/7) can be written as 



y(r) 



2 Cf 



dq 



sm qr 
qr 



av{q) + 



' , / singr 

dq 

V 9?" 



aviq) 



(4.1) 



'0 \ -rf / Jflf 

At short distances, we can expand the singr factor in the low momentum integral, which 
gives the IR power corrections. Making the further assumption that ayiq) has no large 
power corrections at large q and may be well approximated by its perturbative part 
av\PT{q) above fif 



av{q) 



a 



V\PT 



(q) 



(4.2) 



assume the nonperturbative Fourier transform eq.( |2.7| ) exists at least in a formal sense, in particular 
that a long distance Blogr + C contribution is not present, as previously observed. 
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(this assumption will be modified below, eq. (|4.7| )), one ends up with the r — > expansion 



V{r) ~ VpT{r,fif ) 



2 Cf 



IT 



aviq) 



dq q^av{q)+0{r^) 



where 



(4.3) 



2 Cf 



IT 



dq 



( sin qr 



av\PT{(l) 



(4.4) 



\ qr 

is the IR subtracted perturbative potential |^. The normalization of the standard 0{r^) 
and 0{r'^) renormalon-related power corrections in eq.(4.3) is thus given^ by low-energy 
moments of Qy(g). Note that the 0{r^) term is actually infinite, as expected from the 
divergent IR behavior of ay{q). In particular in quenched QCD eq.(|3.5|) implies for — > 



Cf aviq) 



2 K vr / , vr 
^+12(^ + ^2 



(4.5) 



But since the 0{r^) term contributes only an overall normalization constant to the poten- 
tial, which in this section is left arbitrary, one can drop it out. On the other hand, the 
O(r^) and higher order r-dependent contributions are finite. In particular, using eq.([4.5|) 
as a rough approximation to ay{q) in the range < q < jif one obtains in quenched QCD 
for r — > (ignoring any constant term) 

V{r) VpTir, /x;) + (^/./ + + ^/.)) + 0{r% (4.6) 

Let us now modify the previously mentioned assumption, in order to deal with the 
possibility that a 0{l/q^) power correction is actually present in av{q)- Such a correction 
has been first suggested in |jll|] as a consequence of new physics related to confinement, 
leading to a 0{r) linear correction to the potential at short distances, of the same size (and 
sign) as the standard long distance correction related to the string tension. It should be 
noted however that a short distance linear piece may have a more conventional (although 
still non perturbative) infrared origin, as indicated by the position of the leading IR renor- 
malon present in ^(g), which also suggests @ the presence of a 0{l/q^) correction. Let 
us thus assume that for — > 00 



av{q) ^ Oiv\PT{q) + 



2Ko 1 



(4.7) 



Cf 

with Kq ^ K m general. To deal with this correction, one can use the general method of 
1 12], or more conveniently, introduce a new coupling av{q) (different in general from the 
one in section 2, see below), related to the original av{q) by 



Cpaviq) = Cpaviq) + 



(4. 



^The corresponding expressions in term of non-local operators can be found in the effective field theory 
framework of [h Oil. 
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such that the redefined coupUng aviq) is essentially given by its perturbative part (which 
coincides with that of ay (9)) at large 



dviq) ^ aviPTig) 
with no substantial power corrections. Thus from eq.(^) 

av{q) 



V{q) 



AttCf- 



(4.9) 



(4.10) 



and, upon taking the Fourier transform 



V{r) = Ko r + 6V{r), 



(4.11) 



where 6V{r) is given by eg. (K^) and ( |2.10D , but with ay(g) now defined by eq. ([4.8D . Note 
that for Kq = K this definition coincides with that of section 2 (assuming C = 0, see the 
comment after eq.( 2.13| )). Thus, introducing a factorization scale jif as in eq.(4.1) we have 



5V{r) 



2 Cf 



vr 



, / sm \ , ^ , / sm gr , , , 
dq oiv{q) + dql av{q) 



(4.12) 



Since ay{q) has no large power corrections, it can be approximated by its perturbative 
part oiy\pF{q) above some scale /i/, and one deduces the short distance expansion 



5V{r) ~ VpT{r,fJ.f) 



2 Cf 



vr 



dqdviq)-— dq q dv{q) + 0{r 
6 ^0 



(4.13) 



From eq.( [4.5D and (4.8) we get for ^ 



Thus, dropping again the (infinite) 0{r^) term, and using eq.( 4.14| ) for q < Hf, we obtain 



6V{r)^VpT{r,iif) + 
hence from eq.( 4.11j ) 



288 



(4.15) 



V{r) ~ VpTir,Hf) + Ko r + 



f 2{K-Ko) 
\ Svr 



1 ,3 , ,.4 \ Jl 



288' 



(4.16) 



which of course agrees with eq.( [4.6D for Kq = 0. The correlation between the coefficient 
of the 0[r) correction (which is independent) and that of the standard OPE 0{r'^) 
correction should be noted. For Kq ^ 0, we get a neat derivation of the well-known 



statement [11| that the appearance of a linear short distance term in V{r) is equivalent 
to the presence of a 0{l/q'^) correction in the standard av{q) coupling. Moreover, for 
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Kq = K, one obtains the straightforward, but interesting, result that the appearance of 
the hnear short distance term is equivalent to the statement that the modified coupling 
aviq) of section 2 (rather then av{q)) has no corrections. 

One might attempt an analysis of the lattice short distance data of |13] based on 
eg. ( 4. 11 ) and ( f4.13D . Vpxir, ^f) could be evaluated from eq.(4^) by solving the known |14| 
3-loop renormalization group equation for av\pTiQ) ^^^^ performing the integral, similar to 
the single dressed gluon 'renormalon integral' (with IR cut-off) in [12[ 15|, while the power 
corrections should be fitted. Unfortunately, one finds that the perturbative expansion of 
the V-scheme coupling beta function is not reliable at values of small enough that the 
low momentum integral in eq.( [4.12| ) can be meaningfully expanded and parametrized in 
term of a few power correction terms, even at the shortest values of r presently measured 
on the lattice. Thus no reliable fit of the power corrections can be performed yet. It should 
be noted that in the present approach standard IR power corrections appear from an 
OPE-like separation of long and short distances in the Fourier transform of the momentum 
space potential, and their presence is mandatory. This is to be contrasted with the result 
of I IS], where no power corrections were needed'^ if the potential is predicted in term of 
the renormalization group equation of the position space effective charge ap associated 
|17] to the force F{r) = ^ = Cp-^-^--^- However, the implicit definition^ of the power 
corrections in the later case is different, and does not make use of a momentum space IR 
cutoff to separate long from short distances. 



5. Conclusion 

We have shown that it is possible to fix in a natural way the 0{A) renormalon ambiguity 
in the pole mass, thus giving a nonperturbative definition of the pole mass in QCD at this 
level of accuracy, which represents a natural nonperturbative extension of the 'potential 
subtracted' mass, in the spirit of the IR finite coupling approach to power corrections. 
This definition is an optimal one, in the sense the prescription is to remove from the 
heavy quark potential contribution to the self-energy those terms and only those one (the 
confining ones contained in Vconf {'>')) which would give a meaningless (infinite) result for 
the pole mass. For instance, one should not remove from dV{r) the 0{l/r) 'Liischer term' 
to include it in Vconf (see eq.( p.llD ) (which, moreover, would make the modified IR finite 
V-scheme coupling av{q) non-asymptotically free!). The applications of the proposed mass 
definition are similar to those of the 'potential subtracted' mass, to which it provides the 
leading power correction, allowing an accurate relation to the standard MS mass, but it 
can be used consistently with non-perturbative extensions of the Coulomb static potential 
(such as implied by phenomenological potential models or the potential determined on the 
lattice). The remaining challenge is to fix the ©(A^/m) ambiguities in the pole mass arising 
from higher order renormalons. 

^The alternative analysis of also finds no room for power corrections. 

*The convergence of the expansion of the ap beta function is only slightly better than that of 
the momentum space av beta function, which presumably makes a quantitative analysis of the power 
corrections difficult also in the scheme of |13|. 
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We have also discussed the OPE like analysis of the short distance potential. The 
magnitude of the standard OPE contributions have been estimated from eq.([4.5|). However, 
the resulting short distance expansion is unreliable at the lowest values of r measured so 
far on the lattice, due to the poor convergence of perturbation theory for the momentum 
space V-scheme coupling beta function. 
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